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Abstract 

A pair of new colored particles produced near the threshold can form a bound state and 
then annihilate into standard model particles. In this paper we show how the bound state 
signals can distinguish between two scenarios that have similar particle content and inter- 
actions but different spins, such as the minimal supersymmetric standard model (MSSM) 
and universal extra dimensions (UED). We find, for example, that bound states of KK 
gluons (KK gluonia) have an order of magnitude larger cross sections than gluinonia and 
they may be detectable as resonances in the bb, tt or 77 channels if the KK gluon is very 
light and sufficiently long-lived. KK gluonia can be distinguished from gluinonia by their 
much larger cross sections and distinct angular distributions. Similarly, KK quarkonia 
can be distinguished from squarkonia by the size of their diphoton cross section and 
by their dilepton signals. Since many of our results are largely determined by gauge 
interactions, they will be useful for many other new physics scenarios as well. 
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1 Introduction 

Studying the new particles that will be produced at the LHC will likely be a non-trivial task. 
In many scenarios the particles of the new sector are charged under a new parity that makes 
the lightest such particle a dark matter candidate. The new particles will then predominantly 
be produced in pairs, giving multiple decay products some of which will be undetectable 
(including at least two dark matter particles). The reconstruction of such events with many 
objects and missing energy in the final states will be complicated and often ambiguous. It 
may therefore be useful to also look at signals arising from the near-threshold formation and 
annihilation of QCD bound states of the new particles. The bound states behave as resonances 
which annihilate primarily into just two particles (or jets) and no missing energy, so the 
analysis and the interpretation of the signal are much more straightforward than in the more 
conventional channels. In particular, detecting the peak in the invariant mass distribution of 
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any of the possible annihilation channels provides a direct measurement of the particle mass. 
The disadvantages of the bound state signals are their relatively small cross sections and the 
fact that the dominant annihilation mode is into dijets (which have large QCD background), 
with only small branching ratios for annihilation into cleaner signals such as 77. However, 
whenever the bound state signals are observable, they can provide an entirely independent 
method for characterizing or verifying the properties of the new particles. 

In this paper we study how the spins of the new particles are reflected in the properties 
of their bound states and the resulting signals at the LHC As an example, we compare the 
bound states of the level- 1 Kaluza-Klein (KK) modes in the universal extra dimensions model 
(UED) with bound states of the superpartners in the minimal supersymmetric extension of 
the standard model (MSSM). In both scenarios, the new particles are charged in the same way 
under the standard model gauge groups, and the masses of the particles are to a large extent 
free parameters. The fundamental difference between the two scenarios is the spin of the new 
particles. The KK modes in UED have the same spins as the standard model particles (which 
are the zero modes), while the spins of the superpartners in MSSM are different. However, the 
collider signatures of the decay products of a pair of level- 1 KK modes of UED and a pair of 
analogous MSSM superpartners are unfortunately very similar and it is often challenging to 
determine the spin [IllSlElSlElElIIlElinillDlIIIlinillB]. It is therefore interesting to study 
whether and when the bound state annihilation signals can give information about the spin. 

While much is known about bound states of the MSSM particles (see [H] and references 
therein), the bound states in UED models have not been explored, with the exception of bound 
states of KK quarks in the context of a lepton collider [T5l IT6] . Here we will study the various 
possible bound states of colored particles of UED in the context of the LHC and compare them 
to the corresponding bound states in the MSSM. We will see as we go along that once the spin 
of the particle is specified many of the results are largely model-independent since they are 
determined by gauge interactions. Thus our comparisons between particles of different spins 
(in the same color representation) will hold even more generally. 

We start in section[2]by discussing the UED model in the context of our study and reviewing 
the general formalism for bound state computations. In sections [3HS] we present the production 
cross sections, branching ratios and angular distributions for the various bound states and 
compare signals obtained from UED and MSSM. In section [6] we simulate the bound state 
signals and the relevant backgrounds in the dijet, bb, ti, 77 and e~^e~ channels and estimate 
the LHC reach for the various cases. We summarize our conclusions in section [71 

2 Setup 

2.1 Particle spectrum of UED 

In the simplest UED scenario [TT] (for a review, see [H]), all the standard model fields are 
propagating in a single extra dimension of size R ~ TeV~^. The right- and left-handed 
standard model fermions are represented by separate fields in 5 dimensions, and each has 
two chiralities when reduced to 4 dimensions. To restrict each fermion zero-mode to a single 
chirality, the extra dimension is assumed to be an j^L^ orbifold. The zero modes of the 
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unwanted chiralities of the fermions and the 5th components of the gauge fields are then 
projected out by declaring them to be odd under the Z2 orbifold symmetry. The KK modes of 
the right- and left-handed quarks will be denoted by g|j and q'l, or collectively by q*. We will 
refer to and as right- and left-handed KK quarks, even though each one of them is a full 
Dirac fermion. (Similarly, we will refer to the MSSM partners of the right- and left-handed 
quarks, as right- and left-handed squarks.) All the standard model interactions are contained 
in the bulk Lagrangian. Since in 5 dimensions the gauge, Yukawa and quartic-Higgs couplings 
have negative mass dimension, this is an effective theory with a UV cutoff A which is assumed 
to be above the compactification scale: A > 1/R. 

In the 4-dimensional description, the theory contains towers of KK modes for each standard 
model particle, with the mass of the n-th mode given (at tree level) by 

ml = ^ + ml (2.1) 

where mo is the mass of the zero mode (that is the standard model particle itself). As a result, 
all the KK modes at a particular level n are approximately degenerate, with masses ordered like 
in the standard model. Loop corrections shift the masses [T9l [20] , and the resulting spectrum 
of the first KK level has the KK gluon as the heaviest particle, followed by the KK quarks, 
and then the KK excitations of the electroweak gauge bosons, the Higgs, and the leptons. The 
lightest KK particle (LKP), which is typically one of the neutral KK gauge bosons, is a dark 
matter candidate [21], [22] . At tree level, the LKP is stable because of momentum conservation 
in the extra dimension. When loops are taken into account it remains stable because the 
theory retains KK parity under which odd-n modes are charged. The KK parity also prevents 
tree-level contributions to the electroweak observables, thus allowing a compactification scale 
1/i? as low as a few hundred GeV. For our purposes it is important that due to KK parity 
the level-1 KK particles will be produced in pairs, and thus they can form bound states. 

In general, the theory also includes boundary terms (on the orbifold fixed points) whose 
coefficients are free parameters. The boundary terms can be assumed to be symmetric under 
the interchange of the two orbifold fixed points, and then the theory still preserves KK parity. 
However, the mass spectrum of the KK modes can be affected dramatically. Even if bound- 
ary terms are absent at a particular scale, they are generated by the renormalization group 
running [23l [T9l [20] . and the resulting corrections to the masses are of the order of the loop 
corrections. The spectra of [191 120] described above were obtained with the simplifying as- 
sumption that the coefficients of the boundary terms vanish at a specific cutoff scale. However, 
as has been shown in [21] for the case of the gauge kinetic term, even when the coefficient of 
the boundary term is not much larger than the expectation from naive dimensional analysis, 
there is a large effect on the spectrum of the KK modes. Similarly, the effects of boundary 
kinetic and mass terms for a massive scalar field have been analyzed in [25], where the results 
were applied to the electroweak sector and it was found that the identity of the LKP was 
sensitive to the boundary terms. 

In both UED and MSSM, annihilation may or may not be the dominant decay mode of the 
various bound states, depending on the intrinsic decay rates of the constituent particles, which 
in turn depend on the mass spectrum of the model. Since the theoretical and experimental 
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constraints on the possible mass spectra are a question that is decoupled from the bound 
state analysis, we will leave it out of the scope of this paper. We will assume that the 
bound states decay predominantly by annihilation, but the reader should remember that if the 
constituent particles are not sufficiently long-lived, the cross sections will need to be multiplied 
by appropriate branching ratios based on the annihilation rates that we compute here and 
the model-dependent single-particle decay rates. Note, however, that annihilation branching 
ratios close to 1 are not implausible. In particular, in the MSSM there exist various motivated 
scenarios in which this happens to be true for the gluinonium [2] or the stoponium Or 
looking at this from a different perspective, the observation or non-observation of bound state 
signals can give us certain information about the spectrum of the model. 

In UED, the spectrum of the level-1 KK modes plays a crucial role in determining which 
of them are sufficiently long-lived for their bound states to decay by annihilation rather than 
by the decays of the constituent particles. If the assumption of [T^ [2D] that the boundary 
terms vanish at the UV cutoff were true, the KK gluon would have strong two-body decays 
into a KK quark and an antiquark, and the KK quarks could decay electroweakly into a KK 
electroweak gauge boson and a quark. In this case, similarly to what happens in much of 
the parameter space of the MSSM [T3j, the branching ratios for annihilation will be small. 
However, we believe (although without constructing explicit examples) that the presence of 
boundary terms can change the spectrum in ways that would make these decays kinematically 
forbidden for some of the particles. 

For example, if the KK gluon g* becomes lighter than the KK quarks g*, its 2-body 
decays will be replaced by 3-body decays into a KK electroweak gauge boson, a quark and 
an antiquark through diagrams involving an off-shell KK quark, a process suppressed by the 
electroweak coupling and the KK quark mass. This is analogous to the decay of the gluino into 
a neutralino, quark and antiquark in MSSM scenarios with squarks heavier than the gluino, 
which easily makes the gluino sufficiently stable for our purposes |T3]. Similarly to the case of 
the gluino, there is no need for an unnaturally large gap between the KK gluon and KK quark 
masses in order for the annihilation to dominate. For example, suppose that the dominant 
decay process is g* — )■ qqB*, where B* is the KK hypercharge gauge boson. Then the rate is 

Tg* ^ — a, — — — 1 [^] rn^* ~ 10 M 1 [^] (2-2) 

where for simplicity we assumed nig* ^ nig* ^ uib*^ On the other hand, the bound state 
annihilation rates are 

Tannih ~ «f rUg* ~ 10"^ ITlg* (2.3) 

which can easily dominate over 2rg* , especially when we include the numerical prefactors that 
appear in the exact expressions for Fanmh (see appendix IA.1|) . which are as large as O (100) 



for color-singlet bound states, primarily due to multiple powers of color factors. 

The collider signatures of the 3-body decays of the KK gluon or gluino have been studied 
in [9] in an attempt to distinguish between UED and MSSM. One of the goals of the present 



^Thc assumption riig* 3> nig* is not essential for the validity of this estimate. We have checked that for 
niqi. ~ 2mg* , for example, the exact result differs from (|2.2p by less than a factor of 2. 
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paper is to find out to what extent detecting the annihilation decays of KK gluonia (which 
are bound states of two KK gluons) compared to gluinonia (bound states of two gluinos) can 
provide an additional method for distinguishing between UED and MSSM and determining 
the properties of the underlying particles. While gluinonia have been studied extensively in 
the literature, our paper is the first study of KK gluonia. 

Another object of our study is KK-quarkonia, which are KK quark-KK antiquark bound 
states. These have been studied in [151 [ISj in the context of their production and detection at 
a lepton collider. It was found that the spectrum of [T^ I2U] does allow the KK quarks to be 
sufficiently stable for forming KK-quarkonia, but not for the annihilation decays to dominate 
over the single KK quark decays|§ Including boundary terms can probably change these results 
in either direction, and we are particularly interested in the situation in which some of the 
KK quarks are more long-lived so that the branching ratio for annihilation is enhanced. 

Other than that, it is useful to note that many of our results are largely determined by 
gauge interactions alone and depend only on the properties of the binding particles rather 
than the full spectrum of the model and are therefore valid much more generally than just for 
MSSM or UED. It is therefore useful to study all the possible bound states in these models, 
even if some of the binding particles are not sufficiently stable in generic MSSM or UED 
scenarios. With this motivation in mind, our study covers all the possible bound states of 
pair-produced colored particles in MSSM and UED. 



2.2 Bound state formalism 

Assuming that the masses of the pair of particles satisfy mi,m2 ^ Aqcd, their dynamics can 
be described by the single-gluon exchange potential with the Coulomb-like form 

V{r) = -Cy (2.4) 

where Og denotes the strong coupling constant evaluated self-consistently at the scale of the 
inverse Bohr radius = Ca^/i, where /i = 17111712/ {mi + m^) (while the plain will be 
reserved for its value at the scale of mi and m^) and 

C = ]^{Ci + C2- qi2)) (2.5) 

where Ci and C2 are the quadratic Casimirs of the color representations of the two particles 
and C{i2) of the bound state. For bound states considered in this paper, the values of 
are between 0.11 and 0.15, so using (12.40 should be a good approximation, even though 
subleading corrections may have sizeable effects and it would be desirable to compute them, 
along with higher-order corrections to the production and annihilation processes, in a future 
work. Table miists all the possible pair production processes in UED and MSSM and specifies 
in each case which of the color representations of the pair have an attractive potential and 
which do not, based on the sign of C. The values of C for all the attractive configurations 

■^It was claimed in [TS] that KK top quarks can be very stable so that their annihilation decays dominate, 
but this was incorrect, as explained in [27) . 
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UED 


MSSM 


binding 


non-binding 


99 -> 9*9* 


99 99 


1, 8s, 8a 


10,10,27 


qq g*g* 


QQ -> 99 


1, 8s, 8a 




gg q^Tt 


99 QQ* 


1 


8 


qq q*q* 


qq — )■ qq* 


1 


8 


qq — )■ q*q* 


qq —J- qq 


3 


6 


qg q^g" 


Q9 Q9 


3,6 


15 



Table 1: Pair production processes and the color representations of the pair. Two color-octets 
can form an octet with either a symmetric (oc d"'^'^, denoted 8s) or antisymmetric (oc /"'"^, 8a) 
wave function. The Clebsch-Gordan coefficients for all the color decompositions can be found 
in [28]. We have not explicitly listed processes obtained by replacing particles by antiparticles. 



UED 


MSSM 


SU(3) 


C 


(9*9*) 


m 


1 


3 






8 


3/2 


(gV) 


{qq*) 


1 


4/3 


(gV), 


{qq),{q*q*) 


3,3 


2/3 


iQ*9*),{¥9*) 


{Q9), {q*9) 


3,3 


3/2 






6,6 


1/2 



Table 2: Bound states in UED and MSSM, their color representations and the strength of 
their potential fl2.4p . The possible spins of these bound states will be discussed later. 



are given in table |21 In cases where the bound states are colored, they will further hadronize 
with ordinary quarks or gluons to become color-neutral (if they are sufficiently long-lived for 
this to happen). However, these processes will be happening at much larger distance scales 
than both the production and annihilation of the bound states and we therefore ignore their 
effects. 

One can get the matrix elements for bound state production and annihilation by represent- 
ing the bound state as a superposition of states of two free particles with momenta distributed 
according to the bound state wavefunction ip{r). For an 5*- wave bound state, neglecting the 
dependence of the short-distance process on the momenta, it can be easily shown that the 
matrix element between the bound state and the particles from which the pair is produced is 
given by [23 ED] 

7Wbound=^-Mo (2.6) 

where A^o is the matrix element describing the production of the free constituent particles at 
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threshold and /i is their reduced mass. The cross section will then be proportional to 

l.(0)P^^. (2.T) 

More specifically, the bound state production cross section can be written in terms of the 
near-threshold production cross section of the pair of particles ao{s) a^ 

c^bound(5) = ^^1^(0) p27r 5(5- M^) (2.8) 

where M ~ mi + m2 is the mass of the bound state and 



is the factor that makes the continuum production cross section ao{s) vanish at threshold. For 
mi = 7712, P is the velocity of each particle in the center-of-mass frame. 
The annihilation rate of the bound state into two mass-mo particles is 



iOW L ml 



647rmim2\/ (M/2) 



d9 sine 



1 for identical \ / 1 for identical \ 
X- X- (2-10) 

2 bound particles I \ 2 final particles / 

where A4o{0) is the matrix element between the pair of particles that form the bound state 
(with a particular spin and color representation) and the annihilation products, for any po- 
larization and color state of the bound state, and the sum is over the colors and polarizations 
of the products. 



3 KK gluonia vs. gluinonia 

In this section we will study bound states of pairs of color octets: KK gluons g* (spin 1) in 
UED and gluinos g (spin 1/2) in MSSM. Later in this section we will also look at bound states 
of spin color octets. 

It is useful to classify the various possible bound states in the model according to their 
spin J, color representation, parity P, charge conjugation C, and the transformation under the 
chiral U{1) symmetry of the quarks. This allows one to immediately determine the possible 
production and annihilation channels and describe the corresponding processes using effective 
interaction vertices involving the bound state and Standard Model fields. These effective 
vertices can then be used to simulate the bound state processes in an event generator. The 
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color 




can couple to 




1,8s 


0++ 


pQ- ^ 




8a 


1+- 






1,8s 


2++ 




{.99) 


1,8s 


0-+ 


GP'^Gp^, li-fH 




8a 


I— 





Table 3: KK gluonia {g*g*) (UED) and gluinonia (gg) (MSSM) and their possible couplings 
to gluons and quarks. Here G^'^ is the gluon field strength and G^'^ = \e^'^'^"Gpa\ q denotes 
any quark fiavor, while t refers to the top quark whose mass we do not neglect. In part of the 
spin-1 cases, the coupling is not to the bound state operator but to its derivative d^Vy. 

coefficients of the vertices can be determined by matching to the relevant Feynman diagrams 
of the free pair of particles (times the wavefunction at the origin and the other factors). 

In table [3] we list the allowed KK gluonia and gluinonia (taking the spin-statistics theorem 
into account) and the effective vertices through which each of them can couple to gluon or 
quark bilinears. The possible couplings to a pair of photons (for color-singlets) are like the 
couplings to a pair of gluons. 

To determine the charge conjugation properties of the color-octet bound states, it is useful 
to know that the action of charge conjugation C on a non-abelian gauge field is [311 |32l 133] 

where are the generators of gauge transformations. The minus sign here can be described by 
saying that gluons have C = — 1 (this is similar to the photon which transforms as — )■ — A^, 
thus having C = — 1 in a more straightforward sense). Analogous transformation rules apply 
to KK gluons and gluinos. In the same sense, we find that C = +1 for 8s KK gluonia and 
gluinonia and C = —1 for 8a KK gluonia and gluinonia0 

From table [3] we see that the scalar KK gluonia and pseudoscalar gluinonia can be produced 
only by gluon fusion (but can annihilate also to tt). The spin-1 KK gluonium cannot be 
produced at all from gluons or massless quarks while the vector gluinonium can be produced 
from quarks. The tensor KK gluonium couples to both gluons and quarks. We have checked 
that the explicit results that we will now present, based on diagrams in figures [T] and 121 indeed 
match these effective vertices. The tensor KK gluonium happens to couple to gluons only 
through GP^Gpy but not GP^D^D^Gp^. 

While ()2.6p - ()2.7p assume the binding particles to be distinct, the relation (|2.8|) . with |V'(0)P given by 
()2.7p . is valid also if they are identical. Later, in (|2.10[) . |?/'(0)p again refers to the expression p.7p even in 
case of identical particles. 

^Our result for the charge conjugation of the 8a gluinonium {{gg), 8a, J^^ ~ 1~~), perhaps just the 
convention for defining the sign of C that would describe (|3.ip , differs from that of [34l |35l |36] . Note that 
our result (but not theirs) allows the coupling of this gluinonium to massless qq pairs (via the vector current 
q^^q), which must be possible as we know explicitly from the diagrams. 
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Figure 1: Diagrams for production or annihilation of a pair of KK gluons. For bound states, 
the diagrams with s channel gluon do not happen to contribute. 




Figure 2: Diagrams for production or annihilation of a pair of gluinos. 



The parton-level bound state production cross sections can be written as 

^hound,ij{s) = Pi, C(3) n^alal 6{s - M') (3.2) 

where ({3) = Yl'^=i ^/n^ — 1-2 takes into account the contributions of the radial excitations 
and the prefactors Pij for producing the various bound states from partons i,j are presented in 
table m for KK gluonium and gluinonium. These results are based on the diagrams in figures [1] 
and m and factors such as 02. 7p . More details are given in appendix |X1 After convoluting fl3.2p 
with the parton distribution functions a^ 

C(3)7r2a3c,2 ri dx^ ..^ [M^ 



'^bound 

S 



-^^ys X ' ' \xs J 



and multiplying by the appropriate branching ratios, we obtain the cross sections for dijet, 
bb and ti final states as shown in figure [3l We also show to what extent the results depend 
on the mass ratios mg*/mg* or rriq/mg by varying them between 1 and oo|§ We see that the 



^We are using NLO MSTW 2008 PDFs [37] evaluated at the scale M/2, and the center-of-mass energy 
= 14 TcV for the LHC. 

^In practice, the ratio must be somewhat larger than 1 in order for the annihilation rates to dominate, and 
the extreme limit of rriq* /nig* oo is unphysical for UED (although it may be relevant for some other theory 
that docs not contain KK quarks). 
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process 


J = 0,J, = 


J=l,\Jz\=l 


J = 2,\J,\=2 


J 


= 2, J,| = 1 


99 (^V)i 
99 {9*9*)8s 


729/128 
729/512 




243/8 
243/32 




99 {99)1 
99 {9 9) 8s 


243/64 
243/256 








99 (00)1 
gg (00)8s 


243/128 
243/512 








(^V)i 

{9*9*)8s 








M 

5 


r 2m2, ' 
( 2m^. y 


QQ ^ {9 9) 8a 




9 /m? - 

4 Vm? + m?-) 







Table 4: KK gluonium {g*g*) (UED), gluinonium (gg) (MSSM), and octetonium {(fxp) produc- 
tion cross section prefactors Pjj of f l3.2p . The different columns correspond to different values 
of the spin J of the bound state and its polarization (the projection on the beam axis). 

signals will typically be an order of magnitude larger for KK gluonia than for gluinonia of the 
same mass (except if the KK quarks are very heavy, in which case the bb and ti signals of KK 
gluonia disappear). 

It is interesting to ask what creates this large difference in the cross sections. Unless the 
KK gluons (or gluinos) are heavier than about a TeV, production via the gg channel dominates 
due to its higher luminosity, so let us discuss it first and understand the order-of-magnitude 
difference seen already in the parton level expressions of table HI The differences should 
be attributed to a large extent to the more numerous spin possibilities for the KK gluonia. 
According to table HI in the gg channel (in both 1 and 8s representations), the parton- level 
production cross section for J = KK gluonium is only a factor of 3/2 larger than that for 
J = gluinonium. However, the KK gluonium can also be produced in J = 2 state, whose 
cross section (summed over the spin projections) is larger than that of the (J = 0) gluinonium 
by a factor of 8 (or larger than that of J = KK gluonium by a factor of 16/3 ~ 5). Overall, 
the gg channel produces 9.5 times more KK gluonium than gluinonium^ Note that these 
results are determined by gauge interactions alone, and therefore apply much more generally 
than to just UED and MSSM. In particular, the production diagrams do not involve any other 

^In this case, this happens to be exactly the ratio of the near-threshold pair production cross sections 
(despite the fact that for the purposes of bound states we multiply this quantity by a different 17/1(0)!^ for each 
of the attractive color representations and exclude the repulsive ones). In general, this does not need to be 
the case. For example, in the qq channel the ratio between bound state and near-threshold pair production 
cross section is 3 times bigger for the KK gluonia than for the gluinonia because the color representations in 
which KK gluon and gluino pairs are produced at threshold are different. 
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Figure 3: Dijet (top left), bb (top right) and ti (bottom) signals from KK gluonium annihilation 
(solid lines) vs. gluinonium annihilation (dashed lines) as a function of the resonance mass 
(M = 2mg* or 2m^) at the 14 TeV LHC The different curves correspond to different ratios 
of mq*/mg* or rrig/mg that are indicated next to them. The bb signal of (gg) vanishes for 
rriq = rrig, and the bb and ti signals of {g*g*) vanish for mgi/mg* — )■ oo. 



new particles besides the KK gluons or gluinos. 

In the qq channel, the comparison between KK gluonia and gluinonia is sensitive to the 
masses of the KK quarks (relative to the KK gluon) or squarks (relative to the gluino). If the 
KK quarks are very heavy the KK gluonium cross section goes to zero, while there is no such 
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Figure 4: Angular distributions in annihilation of KK gluonia and gluinonia into dijets (top 
left), bb (top right) and ti (bottom) at the LHC for M = 800 GeV (i.e., nig. = m-g = 400 GeV). 
Here 6 is the angle between the beam axis and the direction of motion of the annihilation 
products in the center-of-mass frame. The different curves correspond to different ratios of 
nig* /nig* or nig /nig that are indicated next to them (except for the bb signal of (gg) where 
there is no such dependence). 



effect for the gluinonium. This is because the gluinonium can be produced via a diagram with 
an s channel gluon that does not involve squarks, while a similar diagram for KK gluonium 
vanishes at threshold. On the other hand, the same s channel diagram of gluinonium interferes 
destructively with diagrams with t and u channel squarks, and this makes the gluinonium cross 
section vanish if the squarks are degenerate with the gluino, while there is no such effect for 
the KK gluonium. Note however that the presence of this effect for the gluinonium depends 
on the existence of squarks in MSSM, and thus will not be a general feature of new physics 
models with color-octet spin- 1/2 particles. Similarly, the production of KK gluonia through 
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Figure 5: Simulated angular distributions in annihilation of KK gluonia and gluinonia into bb 
at the LHC. 



the qq channel depends on the existence of KK quarks, and thus may not be present in some 
other model with color-octet spin-1 particles. Another difference between gluinonia and KK 
gluonia is that even though in both cases the bound states are produced with spin component 
±1 along the beam axis, it is an 8a {J = 1) state for gluinonium and 1 or 8s {J = 2) state 
for KK gluonium. The production of the 1 state (with its large color factor (Ci/Cs)^ = 8) in 
UED makes the cross section larger. 

We also see from table H] that KK gluonia will be produced predominantly in spin-2 states. 
If the angular distributions of the annihilation products can be measured, they can be used 
to distinguish the KK gluonium from the gluinonium (which is produced predominantly in 
spin-0 states, and also in the spin-1 state which is important for some of the annihilation 
channels). The angular distributions in the dijet, bb and ti channels are plotted in figure HI 
The curves measured in experiment will be further affected by detector acceptances, effects 
of QCD radiation, mistakes in reconstruction, cuts, and uncertainties in modeling the QCD 
background which needs to be subtracted from the data. However, the remarkable differences 
between KK gluonium and gluinonium that we see in figure H] will hopefully remain. 

To simulate the effects of some of the experimental factors on the angular distributions, 
we generated events with Pythia, and the result (for example, for the bb channel) is shown in 
figure El The details of the simulation and the imposed cuts are the same as will be described in 
section |6] (except that the cut on cos^ is not included and the pr cut is relaxed to pr > 0.3M). 

Another channel we can look at is the annihilation of color-singlet KK gluonia or gluinonia 
into a pair of photons. Even though these processes can only proceed through loop diagrams 
because the KK gluons and gluinos are neutral, it is still interesting to consider this signal 
because the background in the diphoton channel is much more favorable than in the dijet, bb 
or ti channels. As long as the squarks are not much heavier than the gluino, the annihilation 
rate of the color-singlet gluinonium into 77 is ~ 10~^ of its annihilation rate into gg (the 
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Figure 6: Diphoton annihilation signal of gluinonium (thick dashed line), KK gluonium (esti- 
mate; thick solid line), squarkonium (thin dashed line) and KK quarkonium (thin solid line) 
at the 14 TeV LHC as a function of the resonance mass. For gluinonium we assumed rrig = rrig 
(and similarly for KK gluonium). For squarkonium and KK quarkonium we present the contri- 
bution of a single flavor and chirality, and we assumed the charges of the constituent particles 
to be \Q\ =2/3 (for \Q\ = 1/3, the cross section is 16 times smaller). 



exact expression from [36] is given in (1A.35P ). The relevant loop diagrams for KK gluonium 
have not been computed, but for our estimates we will assume that the diphoton branching 
ratio for spin-0 and spin-2 color-singlet KK gluonia is the same as for gluinonium. We expect 
this to be correct up to an (9 (1) factor since the relevant coupling constants and the gauge 
quantum numbers of the relevant particles are the same in both cases. The signal cross 
sections are shown in figure [61 The angular distributions in this channel may also be useful 
for discrimination since for gluinonia only the J = state contributes, while for KK gluonia 
the contribution may be coming from both J = and J = 2. 

We can also compare the KK gluonia and gluinonia (which are bound states of spin-1 and 
spin-1/2 color octets) with bound states of scalar color-octets. Pairs of scalar color-octets 
will be produced by gluon fusion via the same diagrams as the KK gluons in figured! although 
in practice only the quartic vertex is relevant to bound states - "octetonia" (dxp) in 1 and 8s 
color representations. Their production cross sections are included in table Hlfl They are twice 
as small as the gluinonium cross sections. The octetonia decay rates are given in appendix lA. 31 
Since the octetonia are scalars they will annihilate isotropically, and unless is charged under 
the electroweak group or interacts with some other particles, the only significant annihilation 
signal of octetonia will be gg dijets. It will be easy to distinguish them from bound states of 
spin-1 color-octets (like KK gluonia) because of the factor of ~ 20 difference in the cross section 
and the very different angular distributions. It will be more difficult to distinguish them from 
bound states of spin-1/2 color-octets (like gluinonia) because the angular distributions are also 

^The color-singlet result can be found in [38], and it agrees with ours. 
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Table 5: KK quarkonia {q*q*) (UED) and squarkonia {qq*) (MSSM) and their possible cou- 
plings to gluons, quarks and leptons. Pr^l = (1 =t 7^)/2- Antiparticle bound states, where 
relevant, are also present even if not listed in the table. In part of the spin-1 cases, the coupling 
is not to the bound state operator but to its derivative d^V^. 

almost isotropic in the spin- 1/2 case and the difference in the cross section is only a factor 
of 2. If the difference in the cross section is to be used, potential multiplicity of degenerate 
color-octets and higher order QCD corrections will be important. Part of the QCD corrections 
for gluinonium have been computed in [39l [36] ■ 

4 KK quarkonia vs. squarkonia 

In this section we study bound states of fundamental-antifundamental pairs of KK quarks 
(spin 1/2) in UED and squarks (spin 0) in MSSM. Later in this section we will also look at 
bound states of spin 1 particles in the fundamental representation. 

Table |5] classifies the KK quarkonia and squarkonia and their possible couplings to gluons 
and quarks (via the strong interactions) and leptons (via the electroweak interactions). The 
possible couplings to a pair of photons are like to a pair of gluons. The diagrams through 
which the various processes can be realized are shown in figures [7] and [HI 

For both KK quarkonia and squarkonia, the dominant products of the gluon fusion channel 
are spin-0 bound states in which the KK quarks or squarks have same flavors and chiralities. 
These are also the only bound states that can annihilate into a pair of photons, which is the 
most promising detection channel as we discuss in the following. The angular distributions in 
the 77 channel will be isotropic in both UED and MSSM but the size of the cross section can 
still be used for discrimination. 

In the qq fusion channel, the most interesting products are spin-1 bound states of KK 
quarks with same flavors and chiralities since they can annihilate into a pair of leptons, while 
there is no such process for squarkonia. A disadvantage of this production mechanism (which 



15 




Figure 7: Diagrams for production or annihilation of a KK quark-KK antiquark pair. For 
bound states, the diagrams with s-channel gluon do not actually contribute. 




Figure 8: Diagrams for production or annihilation of a squark-antisquark pair. For bound 
states, only the first diagram in each row actually contributes. 

goes through the diagram with a t channel KK gluon from figure [7]) is that it only has access to 
KK quarks of the flavors that are present in the colliding protons, so it will not be useful in the 
likely case that the lightest KK quark is a KK top. This leads us to consider also subleading 
production mechanisms which do not suffer from this problem. In particular, while spin-1 
KK quarkonia cannot couple to a pair of gluons, they can be produced, similarly to J/ tp and 
T [lO], in the process 

99 9{q*x^)i,j=i , x = Lor R (4.1) 

which is flavor-universal. 

16 



process 



J = , = 



J = 1 , J - 



99 {QlQl) 
or (qy^) 



99 {(IlqI) 

or 



or (g^gj^'^) 
orjg^^'^) 



2 

81 



128 
2187 



2m^ 



2 + ^ 

a 



64 
2187 



2m^ 



64 / 2m^* 
2187 V m\+m^* 



512 



2mgmq 



2187 Vm? + m? 



Table 6: KK quarkonium {q*q*) (UED) and squarkonium (qq*) (MSSM) production cross 
section prefactors Pij of fl3.2p . All the bound states are color-singlets. We use q and q' to 
refer to two different flavors of quarks (and similarly for KK quarks and squarks), while the 
notation q^'^ means that the flavor can be either the same or different from that of q. All 
the numbers refer to a single choice of flavors and chiralities. For simplicity, we assumed the 
masses of the different flavors and chiralities of the KK quarks or squarks to be the same. 

The gg and 77 KK quarkonium diagrams in figure [7] describe the same processes as one 
would have for heavy quarkonia (see [H] and references therein for the discussion of toponium 
and [m II2] for quarkonia of new heavy quarks). The production mechanism (14. ip is also 
model- independent (it is determined by interactions with gluons). On the other hand, the qq 
channel is dominated by a diagram involving the KK gluon which does not exist for quarko- 
nia. Another important difference is that the dominant annihilation process of spin-1 heavy 
quarkonia will be into pairs of longitudinal W bosons because of the quarks' large coupling to 
the Higgs which is responsible for their large masses |¥T] . 

The production cross sections for KK quarkonia and squarkonia are given by the expres- 
sions in appendices IA.4I and IA.5I and shown in table |6l Overall, the cross sections are about 
2 orders of magnitude below those of gluinonia of the same mass. This is primarily due to 
the smaller color factors both in |?/'(0)|^ oc and in the short- distance matrix element. As 
a result, the dijet and ti annihilation channels, whose cross sections are included in figures |9] 
and [To], are not promising. 

However, the possibility of tree-level annihilation into 77 (which has low standard model 
background) is much more attractive despite the (a^/a^)-suppressed branching ratio of this 
mode. The cross section for the 77 signal is shown in figure [6] as a function of the bound state 
mass. The angular distributions in this channel are isotropic for both KK quarkonium and 
squarkonium but the signal will be twice larger in UED (when comparing equal KK quark and 
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Figure 9: Dijet signal of the various bound states at the 14 TeV LHC as a function of the 
bound state mass (M = 2m), where all the particles are assumed to have the same mass m. 
The contributions of antiparticle bound states are included wherever relevant. On the left, we 
assume that the signal comes from particular flavors and chiralities of KK quarks or squarks 
and present several examples, while on the right we sum over all the possible combinations. 
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Figure 10: ti signal of the various bound states at the 14 TeV LHC as a function of the bound 
state mass (M = 2m), where all the particles are assumed to have the same mass m. For KK 
quarkonia and squarkonia, on the left we present examples of the contribution from KK tops 
and stops of a single chirality (the results for the other chirality are the same), while on the 
right we assume that both chiralities are close in mass and contribute. 



squark masses) due to the twice larger production cross section. This property can be used 
as a discriminator, although it may be important to take higher-order QCD corrections into 
account. Such corrections for the squarkonium (stoponium) have been studied in [l3llll]. Also, 
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Figure 11: Dilepton annihilation signal of KK quarkonium at the 14 TeV LHC as a function of 
the resonance mass. The cross sections apply for any single flavor of leptons. Left: production 
via qq — >■ q*q* (result independent of nig* unless nig* ^ nig*). Right: production via gg — )■ 
q*q*g for m^* = 2niq* (for nig* = niq* the branching ratios would be an order of magnitude 
smaller) . 
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Table 7: Tripletonium production cross section prefactors Pij of (13.21) . 



we have assumed the total annihilation rate of squarkonium to be dominated by gg, while in 
some cases annihilation into pairs of W, Z or Higgs bosons cannot be neglected ^51 HHl H?! [26] . 

The cross sections for the dilepton annihilation channel of spin-1 KK quarkonia are shown in 
figure [TTl The right plot refers to KK quarkonia produced through the subleading process (14. ip 
which has 

^Ml«3«3,m (4 2) 

243m2 [m^J 
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where as in HO 
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7-1 (7-l)V 7(7 + 1)' (7 + 1)^ 



(4.3) 



On the other hand, squarkonia do not give rise to a dilepton signal. 

While squarks and KK quarks are examples of spin-0 and spin-1/2 particles in the 3 
representation, it may be interesting to consider also bound states of spin-1 particles in the 
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Figure 12: Diagrams for production or annihilation of a pair of KK quarks. The second 
diagram is relevant only if the flavors are equal. 



Figure 13: Diagrams for production or annihilation of a pair of squarks. The second diagram 
is relevant only if the flavors are equal. 

same representation. Such vector color-triplets couple to gluons via 



where W^^ = D^W^-D^W^ with = d^-igsA^^T" and G"^^ is the gluon field strength. The 
second term here is chosen such that tree-level unitarity in the production of vector pairs from 
gg is preserved at high energies, like in the situation when is a gauge boson of an extended 
gauge group |1H]. The W^^ particles may also couple to photons, which can be described by 
including the —ieQA^ term in and adding an electromagnetic term analogous to the second 
term above. The diagrams coupling these vector particles to gluons or photons are the same as 
for squark-antisquark pairs, figure [H The s-channel diagram does not actually contribute for 
the bound states. We do not consider production from qq since it would depend on the (model- 
dependent) couplings of the quarks to the vector bosons. We present the cross sections for the 
resulting bound states, "tripletonia" , in table [71 The production cross section from gluons is 
an order of magnitude bigger than that of squarkonia and KK quarkonia. Furthermore, most 
of the tripletonia are produced with J = 2 rather than J = 0, which leads to different angular 
distributions of the annihilation products. More details are given in appendix IA.6I 

5 Di-KK quarks vs. di-squarks and KK quark-KK gluon 
vs. squark-gluino bound states 

The possible diagrams for the remaining bound states are shown in figures [121 HSl [HI [El 
The resulting cross section prefactors are given in tables [HI and [9l For all of these bound states, 
the annihilation will be almost entirely into dijets, without any cleaner channels (even with 
a small branching ratio) to consider. Furthermore, the dijet signal (see figure [9]) is typically 
even smaller than that of the gluinonium (whose signal we will analyze in more detail in the 
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Table 8: Di-KK quark (g*g*) (UED) and di-squark (gg) (MSSM) production cross section 
prefactors Pij of (13. 2p . All the bound states are in the 3 representation. All the numbers refer 
to a single flavor of that mass, and a single choice of the chiralities, and our notation assumes 
that q' 7^ g. There also exist similar processes in which all the particles are replaced by their 
antiparticles. For simplicity, we assumed the masses of the different flavors and chiralities of 
the KK quarks or squarks to be the same. 
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Figure 14: Diagrams for production or annihilation of a KK quark-KK gluon pair. 




Figure 15: Diagrams for production or annihilation of a squark-gluino pair. For bound states, 
the last diagram does not actually contribute. 



next section), and therefore cannot be seen on top of the QCD background in most scenarios. 
The only exception is if the mass spectrum of the model is very degenerate to the extent that 
signals from many different bound states merge into one. This situation is exemplified in the 
right plot of figure [9] which sums over the flavors and chiralities of the KK quarks or squarks, 
which corresponds to the overly optimistic scenario in which all the flavors and chiralities are 
sufficiently long-lived and close in mass within ~ 5% so that the dijet signals of all their bound 
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Table 9: KK quark-KK gluon {q*g*) (UED) and squark-gluino {qg) (MSSM) bound states 
production cross section prefactors Pij of fl3.2p . All the numbers refer to a single flavor and 
chirality. There also exist similar processes in which all the particles are replaced by their 
antiparticles. 

states merge into a single peak in the invariant mass distribution. Similarly, we may further 
sum the curves corresponding to different types of bound states. Even then, the dijet signal 
will be extremely challenging. 

6 Detection prospects 

We will now analyze to what extent the bound state signals discussed in the previous sec- 
tions will be detectable at the LHCj^ To that end, it is informative to consider the exist- 
ing experimental constraints on the masses of the various UED and MSSM particles. Di- 
rect limits from the Tevatron Run I constrain the size of the extra dimension in UED as 
1/R > 300 GeV [SU [55] (the KK modes masses are nii ~ 1/i?). A much stronger limit can 
probably be obtained from the analysis of data available today. The most stringent indirect 
limit, 1/R > 600 GeV, arises from the inclusive radiative B — )■ Xg'j decay [56]. With the 
assumption that the Higgs is not much heavier than 115 GeV, the same limit is obtained 
from electroweak observables based on LEP data [57] . For the MSSM, the gluino mass is con- 
strained by collider searches to nig > 400 GeV, and squarks are constrained to nig > 600 GeV, 
except for the sbottom (which can be lighter than 300 GeV) and the stop (which can even 
be under 200 GeV) [58l|59l|60l[6Tl[62l[63ll64]. It should be remembered though that MSSM 
and UED have multiple free parameters (the soft SUSY-breaking parameters in MSSM and 
the boundary terms in UED), and experimental bounds usually depend on certain arbitrary 
assumptions about them and thus apply in only part of the parameter space (for an example, 
see [65]). We therefore find it useful to consider also particles that are lighter than the bounds 
quoted above. Note also that many of our results do not depend on the full particle content of 

^LHC signals of gluinonium have been also studied in [49l [50l EU EH [M] and squarkonium in [53l Sll EH |44] . 
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Figure 16: Dijet channel: signal-to-background ratio (left) and the luminosity required for Scr 
significance (right) for gluinonium (dashed lines) and KK gluonium (solid line) at the 14 TeV 
LHC as a function of the resonance mass. For KK gluonium we assume m^* = mg* while for 
gluinonium we present two cases. 



MSSM or UED and can be relevant to bound states in other models that are less constrained. 

To simulate the bound state signals and the dominant standard model backgrounds we used 
Pythia (version 8.120) [661 167] with NLO MSTW 2008 PDFs [37] and SISCONE jet algorithm 
(version 2.0.1) [68] with cone radius R = 1 (except for the tt analysis where we used R = 0.5), 
overlap parameter / = 0.75, no pr threshold on stable cones, and an infinite number of passes. 
We selected several particles defined in Pythia and modified their coupling constants and 
branching ratios such that they would behave according to the bound state effective vertices 
from tables |3] and [51 For the purpose of simulation, we pretended all the bound states to be 
color singlets. We used the BSM Higgs for simulating spin-0 bound states, the Z' for spin-1 
bound states (and T for simulating f l4.ip ). and the KK graviton for spin-2 bound states. We 
have simulated the backgrounds without any i^-factors since our signals do not include higher- 
order QCD corrections either. Such corrections to the pair production processes, the bound 
state wavefunctions and the annihilation processes can be large and sometimes even change 
the cross section by a factor of ~ 2. Part of these corrections have already been computed for 
some of the MSSM bound states [391 ESI SSI HI] but none for UED. Our results will need to 
be re-examined once these corrections are known. 

In the analysis of the dijet channel, we require the two hardest jets within |?7| < 2.5 to 
have pt > 2M/5, the scattering angle in the partonic collision frame to satisfy |cos^| < 0.5, 
and the dijet invariant mass to be within ±15% from the mass of the resonance. In the bb 
channel, we require two tagged jets and assume 60% tagging efficiency for 6-jets, 1% mistag 
rate for gluon and light quark jets and 15% mistag rate for c-jets. Based on [03 EO], we believe 
this level of b tagging performance will be realistic at least for resonances with M < 1 TeV. 
This leaves the standard model bb production as the dominant background. Besides the 
tagging efficiency factor, we use the same cuts as in the dijet channel. In the ti channel, 
we consider the standard model top production processes to be the background and use the 
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Figure 17: bb channel: signal-to-background ratio (left) and the luminosity required for 3o" 
significance (right) for gluinonium (dashed lines) and KK gluonium (solid line) at the 14 TeV 
LHC as a function of the resonance mass. For KK gluonium we assume irtq* = rrig* while for 
gluinonium we present two cases. 
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Figure 18: Semileptonic ti channel: signal-to-background ratio (left) and the luminosity re- 
quired for 3(7 significance (right) for gluinonium (dashed lines) and KK gluonium (solid line) at 
the 14 TeV LHC as a function of the resonance mass. For KK gluonium we assume rriq* = nig* 
while for gluinonium we present two cases. 



procedure described in appendix C of [H] to reconstruct the 4-momenta of the two tops from 
their semileptonic decay products (which is the situation where one W from t — Wb decays 
leptonically and the other hadronically) . We then count events in which the ti invariant mass 
is within ±15% from the mass of the resonance, the tops have pt > M/3, and |cos6'| < 0.5. 

The results for the dijet, bb and ti channels, respectively, are presented figuresHH [TTlandfTSl 
The dijet channel is unlikely to be realistic because of the very small signal-to-background 
ratio S/B (which implies high sensitivity to systematic errors in modeling the background). 
In the bb and ti channels the situation regarding S/B is more promising, especially for KK 
gluonium. As for the fluctuations in the background, the luminosity required for the statistical 
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Figure 19: Diphoton channel: signal-to-background ratio (left) and the luminosity required 
for 3(T significance (right) for gluinonium (thick dashed line), KK gluonium (estimate; thick 
solid line), squarkonium (thin dashed line) and KK quarkonium (thin solid line) at the 14 
TeV LHC as a function of the resonance mass. For gluinonium we assumed mq = rrig (and 
similarly for KK gluonium). 



significance of the signal is achievable for rrig < 450 GeV (if rriq ^ rrig) or mg* < 800 GeV 
in the bb channel, and for rrig < 300 GeV (if rriq ~ rrig) or rrig* < 500 GeV in the ti channel. 

In the diphoton channel we consider qq — > 77 and gg — ?■ 77 to be the dominant back- 
groundsill] We select events in which the two hardest photons within |?7| < 1.5 have pt > 
50 GeV and look in the invariant mass window of ±2% around the mass of the resonance. The 
results are presented in figure fT9l Since gluinonium is unlikely to be lighter than 600 GeV (cor- 
responding to a 300 GeV gluino), the luminosity required to see its 77 signal is too high. On 
the other hand, for KK gluonium the signal is promising up to KK gluon masses of ~ 500 GeV. 
For squarkonium, the signal is viable for squark masses < 350 GeV. While typical squarks 
are expected to be at least as heavy as ~ 600 GeV, a stop can be much lighter. Luckily, 
the stop is also the squark that has the largest chance for being sufficiently long-lived so that 
stoponium will decay primarily by annihilation. Several models in which this happens were 
discussed in [26]. The stoponium 77 signal can thus be within the reach of the LHC. Similarly, 
a KK quark with mass < 400 GeV can have a viable 77 signal from KK quarkonium. Here 
we assumed that the KK quarks or squarks that form the bound state have charge Q = 2/3. 
For KK quarks or squarks with charge Q = — 1/3 the cross sections will be 16 times smaller. 

In the dilepton channel, we focus on e'^e~ processes since they will have the best mass 
resolution. We include Drell-Yan processes as the background. We select events which include 
two electrons with |?7| < 2.5 and pt > 60 GeV and look in the invariant mass window of 
±2% around the mass of the resonance. The resulting reach is shown in figure [20] for several 

^°We compute the significance as S/^/B, and present the value of integrated luminosity for which S/ s/B = 3. 

^^There will also be a contribution from 7+jet and dijet events in which jets are misidentified as photons. 
However, for the heavy resonances considered here, this background can be made subdominant by the tight 
photon identification criteria with only a mild reduction of the signal (see |71j for a recent study by CMS). 
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Figure 20: Dielectron channel: signal-to-background ratio (left) and the luminosity required 
for 3(7 significance (right) for several examples of spin-1 KK quarkonia at the 14 TeV LHC as 
a function of the resonance mass. The solid curves correspond to the left plot of figure [11] and 
the dashed curves to the right plot. 



cases of KK quarkonium. The signal is detectable for some KK quarkonia with KK quarks 
corresponding to the light quark flavors if m^* < 500 GeV. On the other hand, for a perhaps 
more likely case that t\ is light, the e^e~ + jet signal of {t\t*^ is much less promising, 
especially since for this case, as in the right plot of figure [TTl we are already considering the 
more favorable situation of mg* = 2mq* (rather than mg* ~ mq*). 

7 Conclusions 

In this paper we have addressed all the possible 5'-wave QCD bound states of pairs of particles 
from the UED scenario and compared them with analogous bound states of MSSM particles. 
KK gluonium and gluinonium were also compared with bound states of adjoint scalars, and 
KK quarkonium and squarkonium with bound states of color-triplet spin-1 particles. 

We have found that if the KK gluon is sufficiently stable (which is likely to be the case if 
it is lighter than the KK quarks), its bound states signals at the LHC will have cross sections 
that are significantly larger than those of the gluinonium of the MSSM, mostly due to the 
larger number of spin states for KK gluonia. Besides the order-of-magnitude difference in the 
size of the cross sections, the angular distributions of the annihilation products can be used 
as a discriminator since KK gluonia will predominantly form with spin 2 (and some spin 0), 
while gluinonia will be mostly spin (and some spin 1). 

The diphoton signal of spin-0 KK quarkonium is twice as large as that of the squarkonium, 
and is potentially detectable if the constituent KK quarks are sufficiently stable and light. This 
difference in the cross section can help distinguish between KK quarks and squarks. Another 
potential discriminator is the dilepton annihilation signal of the vector KK quarkonium which 
might be measurable. 

Despite the fact that bound state cross sections are very small, we found that in some cases 
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the detection of the signal is reahstic. For example, a 600 GeV KK gluon can give rise to a 
1.2 TeV resonance in the bb channel, with S/B 3% and 3a significance reached at ~ 30 fb~^. 
With even more luminosity, the angular distribution can hopefully be extracted and indicate 
the spin-2 nature of the bound state, confirming that the KK gluon has spin 1. Another 
example is a 400 GeV KK top whose bound state will appear as an 800 GeV resonance in the 77 
channel with S/B 10% and 3a significance at ~ 600 fb~^. Furthermore, many of our results 
remain valid for other models that contain pair-produced particles with color representations 
and spins that our study has covered, but which can otherwise be slightly or completely 
different from the standard MSSM and UED scenarios and have weaker experimental bounds 
on the masses, which opens up additional possibilities. 

Overall, our study has demonstrated how processes of bound state formation and annihi- 
lation, which are easy to compute for a given model and easy to reconstruct at the collider, 
provide additional channels for studying new particles at the LHC. 
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A Cross sections, annihilation rates, angular distribu- 
tions 

In this appendix we give the expressions for the parton-level near-threshold pair-production 
cross section o"o(s) which enters the calculation of the bound state production cross section 
(12. Sp . For UED processes, we used Feynman rules from [721 |5l]. We separate the results by the 
color representation of the pair, its spin J, and the spin component Jz along the direction of 
motion of the first parton (where by, e.g., \Jz\ = 1, we will refer to the sum of the contributions 
with = 1 and = — 1 in cases when the two contributions are equal). We consider all the 
attractive color states and all the possible spins that can be produced from incoming quarks 
and/or gluons. We use q and q' to refer to two different flavors of quarks (and similarly for 
KK quarks and squarks), while the notation q^'^ means that the flavor can be either the same 
or different from q. We use the subscript x {^^ Qx) denote the chirality of the particle 
{x = L 01 R) while 3^ refers to the other chirality. Our expressions for the cross sections refer 
to a single flavor and chirality choice for each of the two produced particles (for squarks and 
KK quarks). Processes that are not listed have vanishing cross sections at threshold. 

We also list the annihilation rates of the bound states into various flnal states and the corre- 
sponding angular distributions (where 9 is the angle between the momenta of the annihilation 
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products and the spin quantization axis, in the center-of-mass frame). 

When a t quark appears as an annihilation product we take its mass into account, but 
otherwise the masses of the quarks are set to zero. For simphcity, we assume that all the KK 
quarks have a common mass rriq* and all the squarks have mass rriq. 



A.l KK gluonia 

The cross sections for KK gluon pair production processes in the near-threshold limit, for the 
attractive color representations, are 



aoi99^9'g^;l,J = 0) = '^^^ (A.l 

Zbbmf 



9 



^0 (99 ^ ; 8s, J = 0) = (A.2 



^0 {99 ^ 9*9' ; 1, J = 2, 1 J,| = 2) = (A.3 



h {99 -> ; 8s, J = 2, I J,| = 2) = (A.4 



, 2m2 \ ^ 27ra2^ 



(,5 ^ ; 8s, J = 2. U.l = 1) = ( ;;^^) 1^ (A.6 



The annihilation rates of the KK gluonia into the various final states are 

729 

r {{9*9*)i,J=o ^99) = — a^alnig* (A.7 

729 

r ((^7V)8s,J=o ^99) = ^ alalrrig. (A.8 

r ((.V)m^o ^ tt) ^ --f - -f j ^^4_^. j (A.9 

^^^' ' ^^^^^-^^ ''^-2564 ['-4) U+4-^? J ^"^"^^^ ^^-'^ 

r (((7V)i,J=2 = ^alalnig. (A. 11 

243 

r ((^V)8s,J=2 ^ = — (A. 12 
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18 / 2?Ti^ 
9 / 2m^ \ ^ 

18 m? \ 2m? \ / 2m2 n ' 



r {{9*9nss,j=2 ^tt) = ^ (^1 



9' 

1 + ^ ^ aW.m„ 



nig* ) \ 6m„* ) \ m% + m^* — mi ' 



(A.16) 

The annihilation rates into qq should be further summed over the different quark flavors. 

The angular distributions in annihilation of KK gluonia (in both 1 and 8s representations) 
into gg are given by 

1 dT 5 

J = 2, J, = ±2: (35 + 28 cos 2^ + cos 4^) (A.17) 

1 smddy 256 

1 dT 5 

J = 2, J, = ±1 : t;—^ = T7I (3 + cos 29) sin^ 9 (A.18) 

1 sm 9 d9 16 

1 dV 15 

J = 2,J, = 0: r-^ = T^sin^^ (A.19) 

1 sm 9 d9 16 

In annihilation into qq, for massless quarks we have 

J = 2, J, = ±2 : t;—^ = (3 + cos 29) sin^ 9 (A.20) 

1 sm6'a6' 16 

1 dT 5 

J = 2, J, = ±1 : y;—^ = TTt (2 + cos 2^ + cos A9) (A.21) 

1 sm9d9 16 

1 15 ^ ^ , 

J = 2,J, = 0: sin' 29 A.22 

r sm d9 16 

and for 

j = 2.j.^±2 1 ''r 5 1 



rsin^rf^ 16 1 + 2m2/3m2, 

X (^3 + ^ + (^1 - cos 2^ j sin^ 9 (A.23) 
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J = 2, X = ±1 : 



rsin^d^ 16 1 + 2m2/3m 



X (^2 + ^ + cos 20 + (^1 - cos 40^ (A.24) 
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J = 2, J, 



: 



1 dT 



15 



rsin^d^ 16 1 + 2m2/3m2, 

n „ 11m? 
X sin^ 29 ' * 



(A.25) 



18m2, 



2m? 

* cos 20 



3m^* 



2m2, 



cos46' 



A. 2 Gluinonia 

The near-threshold production cross sections of attractive gluino pairs are 



o-Q {gg ^ gg; 1, J = 0) 
o"o ^ gg; 8s, J = o) 

-^0 (gg gg ; ^A, J = i, l^^l = i) 

The annihilation rates of the gluinonia are 

243 



128m? 



m? + m? 



3m? 



r((^^)i,j=o gg 
r {{gg)8s,j=o gg 



r ((^^) 



SsJ=0 



r ((^^) 



8a,J=1 



tt 



—alalm~g 

243 _o 2 
128 " " ^ 



m? , 
9^W1 
m^ 



mt 



2m? 



m? \m? + m? — m^ 



45 m? 



128 m? 



2m? 



mt 



m? \ m? + m? 



mt 
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64 V m| + m| 
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2m? 



^ _ m£ 
m? 



m^ + m^ — mf 



(A.26) 
(A.27) 

(A.28) 



2 2 2 \ 2 



(A.29) 
(A.30) 

(A.31) 
(A.32) 
(A.33) 
(A.34) 



The branching ratio into 77 via loop diagrams is [36] 
r((^^)i,j=o ->77) 50^2 



Tm)i.J=o^gg) 8l7r2 
where in the last expression we substituted rriq 



m"' \ ( m? 



3 

m? 
= m 



257r2 
"648 



-OL 



(A.35) 



9- 
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The J = 1 8a gluinonium has the following angular distributions in the annihilation into 
qq for massless quarks: 

1 (JX^ 3 

J = 1,J, = ±1: (1 + cos^^) (A.36) 

J = 1,J, = 0: ^-4^ = isiii'^ (A.37) 
i sm d ad 4 

while for tt: 

J = 1,J, = ±1: _^ = ^^ ^cos^^ (A.38) 

TsmOde 8 1 + /2m? V 1 + /m? J ^ ' 

J = 1, J. = : = 7 1 fsi^' ^ + 4 cos^ ^ (A.39) 

Tsm 6'd6' 4 1 + m^/2m?V rnj J ^ ^ 

A. 3 Octetonia 

The near-threshold production cross sections of attractive pairs of scalar color-octets are 

*„ (99 ^ 8s) = ^ (A.41) 
The annihilation rates of the octetonia are 

243 

r ((00)1 ^ gg) = —a^alm^ (A.42) 
243 

r ((00)8s ^ = —ala^,m^ (A.43) 

A. 4 KK quarkonia 

The near-threshold production cross sections of attractive KK quark-KK antiquark pairs are 

^0 {99 ^ ; 1, J = 0) = ^ (A.44) 

*o f 95"> ^ ; 1. J = 1 . J, = O) = ( 74 V ^ ( A.46) 



-^0 (gg^'^ ^ qWl'^ ; 1, J = 1, J, = -l) = (To (gg(') ^ g^^^'^ ; 1, J = 1, J, = 1 
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(A.47) 



The annihilation rates of the KK quarkonia are 



g4 

r j=o ^ gg) = ^ (A-48) 
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?!^V \ (^2 + ^-^) aVm.(A49) 

r ((gx9^)i,J=o ^ 77) = y a^g'«'mg. (A.53) 
The electroweak rates into fermions / (leptons or quarks) are 

r a,) ^ ( (MM) . (A.4) 

where x ^i-nd 77 describe the chirahty of the KK quark and the fermion, respectively, = 1 for 
leptons and 3 for quarks, Q is the electric charge and is the weak isospin, and we assumed 
m| <^ (2mg*)^. More specifically, for annihilation into charged leptons 

16 (x^ 

r j=i -> t^q) = ^^'c^rj^^^m,. (A.55) 

with 

crr = , crl = ^ , Cifi = (g - T^) % = ^ (g - T3 + t3 cot^ 9wY (A.56) 

where one should substitute Q = 2/3, = 1/2 for up-type KK quarks and Q = —1/3, 
= —1/2 for down- type KK quarks. 
The angular distributions for iqLqi}'^)i,j=i — )■ Q'^*-'-* are 

J=1,J^ = 1: ^ = -sin^- (A.57) 

Tsin^d^ 2 2 ^ ' 
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1 dT 3 ,e , , , 

while for j=i — >■ Q'^'^'-' the same expressions hold for opposite signs of J^. For both 

{qlqt'\j=i -> <1Q^'^ and (gjj^^'^)i, j=i ggW, 

J=1,J, = ±1: ^-4^ = |sin2^ (A.60) 

r sm OdO 4 

1 c^r 3 

J=1,J^ = 0: - = -cos2^ (A.61) 



For (g;;g*)i,j=i ^ 



r sin ede 2 



1 rfr 3 , 

1 smadd 4 

while for j=i — )■ ^^^^ the same expressions hold for opposite signs of Jz- 

A. 5 Squarkonia 

The near-threshold production cross sections of attractive squark-antisquark pairs are 

*o (99 -* ?x5; 1.^ = 0) = ]^ (A.65) 



9 ' 3/ 1 



The squarkonium annihilation rates are 
r {{,qxQ*x\j=o 99) = ^oiWs^q (A.67) 



r (fat;),..o ^ ^0 - ^4 (1 - 4) ^1 - 4 ^ 4_ J (A.68) 

2 



16 

r (feC)i.^=o ^ 77) = -alQWm^ (A.70) 
where Q is the electric charge of the squark. 
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A. 6 Tripletonia 

The near-threshold production cross sections of attractive pairs of vector color-triplets are 



ao {gg ^ WW* ; 1, J = 2, | J,| = 2) = (A.71) 



w 

2/ 



ao {gg WW* ■ 1, J = 0) = ^ (A.72) 



The annihilation rates of the tripletonia into gg are 

512 



r mW*),,j=2 ^99) = -r7^ aiaimw (A.73) 

4U5 

32 

r {{WW*)i,j=o ^99) = ^ alalmw (A.74) 

If the vector bosons have electric charge Q they can annihilate into 77 with rates which can 
be obtained by the replacement 

^Q^a^ (A.75) 
The angular distributions of the annihilation products (both gg and 77) are 

1 dV 



J = 2, J, = ±2 
J = 2,J, = ±l 
J = 2,J, = 



r sin e d9 256 
1 dV 



r sin 9 dd 16 
1 dV 



r sin d dd 16 



-^(35 + 28 cos 2^ + cos 4^) 
256 


(A.76) 


(3 + cos 26) sin^ 9 

16 


(A.77) 


15 . 4 ^ 
— sm 
16 


(A.78) 



A. 7 Di-KK quarks 

The near-threshold production cross sections of attractive KK quark pairs are 

/ 2m2 \^ \2 vra^o 

(«, ; 3. J = 1. J, = 0) = (;;^) (;;^ j ^ (A.79) 



/ 2m2 \2 / m2 \ 2 7ra2/5 



















4 + 






2 + 


m\. 


m]. 
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•^0 [qq' qWu ; 3, J = 1, J, = -l) = <To {qq' ^ qnql' ; 3, J = 1, J, = l) 



" ^ '2 + ^ -^p^ (A.83) 



There are also processes in which the quarks and KK quarks are replaced by their antiparticles, 
which have the same cross sections for opposite signs of Jz- 

The bound states annihilate into the same quark flavors from which they were produced. 
The rates are 

2 \ 2 



r -> ^ 4 (;^) (:=f ) (A,84) 

/ \ 9 / 2m?' \ ^ / \ ^ 



The angular distributions for (g* g* 3 — )■ gg*^') are 

1 dV 3 



J=1,J, = ±1: __:^ = rsin2^ (A.87) 

J=1,J. = 0: i-^ = |cos2^ (A.88) 

i smddu 2 



For (g*g^)3 qq, 



For {q'lq}i')s j=i — ^ (and for {q^q}^)^ j=i — gg' for opposite signs of Jz) 

I dV 3 ,9 , 

i smadd 4 

For the corresponding anti di-KK quarks, the angular distributions are the same for opposite 
signs of Jz- 
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A. 8 Di-squarks 

The near-threshold production cross section for attractive squark pairs is 

ao(«'^^.C;3.J = 0) = (^j ^ (A.94) 

There are also processes in which the quarks and squarks are replaced by their antiparticles, 
which have the same cross section. 

The annihilation rate of di-squarks is 

A. 9 KK quark-KK gluon bound states 

The near-threshold production cross sections for attractive KK quark-KK gluon pairs are 



^r9 ; 3, J = -, J. = -- = " \ ; nalP (A.97) 



/ _ 1 l\ 3m * 

(« - i «■ = 2. = -2 j = 32^g.K! 

Processes with q\ instead of gjj have the same cross sections for opposite signs of Jz- There 
are also processes in which the quark and KK quark are replaced by their antiparticles, which 
have the same cross sections for opposite signs of Jz- 
The annihilation rates are 

r (gV)3, j=3/2 ^ qg) = « « V , ^ (^-^^o) 

64 (jng*+mq.) 

r (gV)3, .=1/2 ^ = ^ « « ^-^^ ^ (A.ioi) 

r ((9V)6, j=3/2 ^ 5^?) = li^Wsm,* (A.102) 

3 

r ((gV)6, j=i/2 -> 9^?) = 7 ^ — ^ (A.103) 

64 +^ 
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The angular distributions for {q\g*) — )■ qg (in both 3 and 6 representations) are 

. 3,3 I dV .9 X 

3 1 1 dT .0.9 . 



3 
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dT 


= 2 cos*^ 
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sin 6 dO 
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+r 
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= 6 cos^ 
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= 6 cos^ 
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~2 ■ 
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smOdO 
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dT 


= 2 sin^ 




~2 ■ 
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sin 9 dO 


2 



3 1 1 dT o9 ,9 . 

. 1.1 1 dr X 

•'=2"'' = +2^ r;ta9d«=='"2 f*-!"*^' 

T 1^1 I dT ^6 . . X 

J=-,Jz = — : ^ = cos2- A.109 

2' 2 Tsin^c/^ 2 ^ ^ 

Processes with instead of g|j have the same angular distributions for opposite signs of J^. 
For the corresponding anti KK quark-KK gluon bound states, the angular distributions are 
again the same for opposite signs of Jz- 

A. 10 Squark-gluino bound states 

The near-threshold production cross sections for attractive squark-gluino pairs are 

1 . 1\ img + 9mgY ^^2, 

i-g (rr 

^0 (qg -^qRg;6, J = J, Jz = — ] = t-^ — : r ^"^^ (A.lll) 

\ 2 2/ 16 nig {nig + nig) 

Processes with qi instead of q^ have the same cross sections for opposite signs of Jz- There 
are also processes in which the quark and squark are replaced by their antiparticles, which 
have the same cross sections for opposite signs of Jz- 
The annihilation rates are 

r [{qg)s,j=i/2 -^qg) = ^ f — ■ (A.112) 

32 [m-g + nig) 

1 Tn~Tfi~ 

r {{q~gh.j=i/2 ^ qg) = — (A.113) 

32 [nig + nig) 

The angular distributions for {qng) — qg (in both 3 and 6 representations) are 

2 i sm b' do' 2 
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ao\qg^qR~g-,^.J=-.J. = --]= ^ na^P (A.llO) 

^ ^/ 288 (mg -I- mg) 



J. = —. 



1 dV 



cos 



(A.115) 



The angular distributions for (^l^) — )■ qg are the same for opposite signs of J^. For the 
corresponding anti squark-gluino bound states, the angular distributions are the same, again 
for opposite signs of J^- 



References 



[1 

[2 

[3] 
[4] 
[5 

[6, 

[7 



[9 
[10 



J. M. Smillie and B. R. Webber, "Distinguishing spins in supersymmetric and universal 



extra dimension models at the Large Hadron Collider," JHEP 10 (2005) 069 



arXiv : hep-ph/0507170^ 



A. Datta, K. Kong, and K. T. Matchev, "Discrimination of supersymmetry and 
universal extra dimensions at hadron colliders," Phys. Rev. D72 (2005) 096006| 



arXiv : hep-ph/0509246f Erratum ibid. D72 (2005) 119901. 



A. J. Barr, "Measuring slepton spin at the LHC," \JHEP 02 (2006) 042 



[arXiv : hep-ph/0511 115 



A. Alves, O. Eboh, and T. Plehn, "It's a gluino," \Phys. Rev. D74 (2006) 095010 
[arXiv : hep-ph/06d5~067 



C. Athanasiou, C. G. Lester, J. M. Smillie, and B. R. Webber, "Distinguishing spins in 



decay chains at the Large Hadron Collider," \JHEP 08 (2006) 055 
arXiv : hep-ph/ 0605286j 



C. Athanasiou, C. G. Lester, J. M. Smillie, and B. R. Webber, "Addendum to 
'Distinguishing spins in decay chains at the Large Hadron Collider'," 



arXiv:hep-ph/0606212 



L.-T. Wang and I. Yavin, "Spin measurements in cascade decays at the LHC," 



JHEP 04 (2007) 032,[arXiv:hep-ph/0605296 



C. Kilic, L.-T. Wang, and I. Yavin, "On the existence of angular correlations in decays 



with heavy matter partners," JHEP 05 (2007) 052 , arXiv : hep-ph/0703085 



C. Csaki, J. Heinonen, and M. Perelstein, "Testing gluino spin with three-body decays. 



JHEP 10 (2007) 107, arXiv: 0707. 0014 [hep-ph] 



M. Burns, K. Kong, K. T. Matchev, and M. Park, "A general method for 
model-independent measurements of particle spins, couplings and mixing angles in 



cascade decays with missing energy at hadron colliders," JHEP 10 (2008) 081 



arXiv : 0808 . 2472 [hep-ph] 



[11] O. Gedalia, S. J. Lee, and G. Perez, "Spin determination via third generation cascade 



decays," Phys. Rev. D80 (2009) 035012 arXiv: 0901 .4438 [hep-ph] 



38 



H.-C. Cheng, Z. Han, I.-W. Kim, and L.-T. Wang, "Missing momentum reconstruction 



and spin measurements at hadron colliders," JHEP 11 (2010) 122 



arXiv: 1008.0405 [hep-ph] | 



G. Moortgat-Pick, K. Rolbiecki, and J. Tattersall, "Early spin determination at the 



LHC?," Phys.Lett. B699 (2011) 158 , |arXiv: 1102.0293 [hep-ph] 



Y. Kats and M. D. Schwartz, "Annihilation decays of bound states at the LHC, 



JHEP 04 (2010) 016[|arXiv:0912.0526 [hep-ph] 



C. D. Carone, J. M. Conroy, M. Sher, and I. Turan, "Universal extra dimensions and 
Kaluza-Klein bound states," \Phys. Rev. D69 (2004) 074018[ |arXiv : hep-ph/0312055 



N. Fabiano and O. Panella, "Threshold production of metastable bound states of 



Kaluza Klein excitations in universal extra dimensions," Phys. Rev. D81 (2010) 115001 
[arXiv: 0804. 3917 [hep-ph] 



T. Appelquist, H.-C. Cheng, and B. A. Dobrescu, "Bounds on universal extra 



dimensions," Phys. Rev. D64 (2001) 035002, ^arXiv : hep-ph/0012100 . 



D. Hooper and S. Profumo, "Dark matter and collider phenomenology of universal extra 
dimensions," \Phys. Rept. 453 (2007) 29, ,arXiv:hep-ph/0701197, 



H.-C. Cheng, K. T. Matchev, and M. Schmaltz, "Bosonic supersymmetry? Getting 



fooled at the CERN LHC," \Phys. Rev. D66 (2002) 056006^ .arXiv : hep-ph/0205314| 
H.-C. Cheng, K. T. Matchev, and M. Schmaltz, "Radiative corrections to Kaluza-Klein 



masses, 



Phys. Rev. D66 (2002) 036005narXiv:hep-ph/0204342 



G. Servant and T. M. P. Tait, "Is the lightest Kaluza-Klein particle a viable dark 



matter candidate?," Nucl. Phys. B650 (2003) 391 [iLrXrv:hep-ph/0206071 



H.-C. Cheng, J. L. Feng, and K. T. Matchev, "Kaluza-Klein dark matter," 



Phys. Rev. Lett. 89 (2002) 211301 |arXiv:hep-ph/0207125 



H. Georgi, A. K. Grant, and G. Hailu, "Brane couplings from bulk loops. 



Phys. Lett. B506 (2001) 207, arXiv : hep-ph/0012379 



M. S. Carena, T. M. P. Tait, and C. E. M. Wagner, "Branes and orbifolds are opaque," 
Acta Phys. Polon. B33 (2002) 2355, arXiv : hep-ph/0207056 , 

T. Flacke, A. Menon, and D. J. Phalen, "Non-minimal universal extra dimensions," 
Phys. Rev. D79 (2009) 056009, ,arXiv: 0811. 1598 [hep-ph] ^ 



S. P. Martin, "Diphoton decays of stoponium at the Large Hadron Collider," 



Phys. Rev. D77 (2008) 075002narXiv: 0801 .0237 [hep-ph] 



39 



[27] 
[28] 

[29] 

[30] 
[31] 

[32] 

[33] 

[34] 

[35] 
[36] 

[37] 



E. De Pree and M. Sher, "Kaluza-Klein mesons in universal extra dimensions," 



Phys. Rev. D72 (2005) 097701, [arXiv:hep-ph/0507313 



M. Beneke, P. Falgari, and C. Schwinn, "Soft radiation in heavy-particle pair 
production: all-order colour structure and two-loop anomalous dimension," 



Nud. Phys. B828 (2010) 69 , [arXiv : 0907 . 1443 [hep^ph] 



M. E. Peskin and D. V. Schroeder, An Introduction to Quantum Field Theory. 
Westview Press, 1995. See Section 5.3. 



V. A. Novikov et al, "Charmonium and gluons," Phys. Kept. 41 (1978) 1 



I. V. Tyutin and B. B. Lokhvitskii, "Charge conjugation of non-Abelian gauge fields," 
Sov. Phys. J. 25 (1982) 346, [Izv. Vys. Uch. Zav. Fiz. 4 (1982) 62]. 



N. V. Smolyakov, "Furry theorem for non-abelian gauge Lagrangians," 



Theor. Math. Phys. 50 (1982) 225,. [Teor. Mat. Fiz. 50 (1982) 344]. 



C. P. Burgess and G. D. Moore, The Standard Model: A Primer. Cambridge University 
Press, 2007. 

J. T. Goldman and H. Haber, "Gluinonium: The hydrogen atom of supersymmetry," 
\Physica 15D (1985) 18l[ 



W.-Y. Keung and A. Khare, "Two-gluino bound states," Phys. Rev. D29 (1984) 2657 



M. R. Kauth, J. H. Kiihn, P. Marquard, and M. Steinhauser, "Gluinonia: Energy levels, 
production and decay," \Nucl. Phys. B831 (2010) 285, ,arXiv: 0910 .2612 [hep-ph] , 



A. D. Martin, W. J. Stirling, R. S. Thorne, and G. Watt, "Parton distributions for the 



LHC," Eur. Phys. J. C63 (2009) 189 arXiv: 0901 . 0002 [hep-ph] 



C. Kim and T. Mehen, "Color octet scalar bound states at the LHC," 



Phys. Rev. D79 (2009) 035011, arXiv: 0812. 0307 [hep-ph] 



[39] 
[40] 
[41] 
[42] 



K. Hagiwara and H. Yokoya, "Bound-state effects on gluino-pair production at hadron 



colhders," JHEP 10 (2009) 049 arXiv: 0909 .3204 [hep-ph] 



V. D. Barger and A. D. Martin, "Quarkonium production at pp colliders," 



Phys. Rev. D31 (1985) 1051 



V. D. Barger et al., "Superheavy quarkonium production and decays: A new Higgs 
signal," Phys. Rev. D35 (1987) 3366 , Erratum ibid. D38 (1988) 1632. 

E. Arik, O. Qakir, S. A. Qetin, and S. Sultansoy, "Fourth generation pseudoscalar 
quarkonium production and observability at hadron colliders," 



Phys. Rev. D66 (2002) 116006] [iLrXiv:hep-ph/0208169 . 



40 



S. p. Martin and J. E. Younkin, "Radiative corrections to stoponium annihilation 



decays," Phys. Rev. D80 (2009) 035026 |arXiv: 0901 .4318 [hep-ph] 



J. E. Younkin and S. P. Martin, "QCD corrections to stoponium production at hadron 



colliders," Phys. Rev. D81 (2010) 055006, [arXiv: 0912. 4813 [hep-ph] 



M. Herrero, A. Mendez, and T. Rizzo, "Production of heavy squarkonium at high 



energy pp coHiders," Phys.Lett. B200 (1988) 205 



V. D. Barger and W.-Y. Keung, "Stoponium decays to Higgs bosons, 
[Phys.Lett. B211 (1988)"355l 



M. Drees and M. M. Nojiri, "Production and decay of scalar top squarkonium bound 
states," ^Phys.Rev. D49 (1994) 4595, .arXiv : hep-ph/9312213, 

G. V. Borisov, Y. F. Pirogov, and K. R. Rudakov, "Pair production of exotic particles 
at ^p>p colliding beams," |Z. Phys. C36 (1987) 217| 

E. Chikovani, V. Kartvelishvili, R. Shanidze, and G. Shaw, "Bound states of two gluinos 



at the Tevatron and CERN LHC," Phys. Rev. D53 (1996) 6653 



arXiv : hep-ph/9602249[ 



K. Gheung and W.-Y. Keung, "Split supersymmetry, stable gluino, and gluinonium," 



Phys. Rev. D71 (2005) 015015, [arXiv:hep-ph/0408335 



E. Bouhova-Thacker, V. Kartvelishvili, and A. Small, "Search for gluino-gluino bound 



states," Nucl. Phys. Proa. Suppl. 133 (2004) 122 



E. Bouhova-Thacker, V. Kartvelishvili, and A. Small, "Search for gluino-gluino bound 



states with ATLAS," Nucl. Phys. Proa. Suppl. 152 (2006) 300 



M. Drees and M. M. Nojiri, "Proposed new signal for scalar top-squark bound-state 



production," Phys. Rev. Lett. 72 (1994) 2324, arXiv:hep-ph/ 9310209. 



G. Macesanu, G. D. McMullen, and S. Nandi, "Gollider implications of universal extra 
dimensions," ^Phys. Rev. D66 (2002) 015009, ,arXiv : hep-ph/0201300, 

G. Lin, "A search for universal extra dimensions in the multi-lepton channel from pp 
collisions at = 1.8 TeV," FERMILAB-THESIS-2005-69. 

U. Haisch and A. Weiler, "Bound on minimal universal extra dimensions from 



B X,7," Phys. Rev. D76 (2007) 034014, arXiv : hep-ph/0703064 



1. Gogoladze and G. Macesanu, "Precision electroweak constraints on universal extra 



dimensions revisited," \Phys. Rev. D74 (2006) 093012| |arXiv : hep-p h/0605207| 

J. L. Feng, J.-F. Grivaz, and J. Nachtman, "Searches for supersymmetry at high-energy 



colliders," Rev. Mod. Phys. 82 (2010) 699, arXiv: 0903. 0046 [hep-ex] 



41 



[59] 

[60] 

[61] 

[62] 
[63] 
[64] 
[65] 
[66] 
[67] 

[68] 

[69] 
[70] 
[71] 
[72] 



CMS Collaboration, "Search for supersymmetry in pp collisions at 7 TeV in events with 



jets and missing transverse energy," Phys.Lett. B698 (2011) 196 



arXiv: 1101 .1628 [hep-ex] | 



ATLAS Collaboration, "Search for supersymmetry using final states with one lepton, 
jets, and missing transverse momentum with the ATLAS detector in ^/s = 7 TeV pp 



collisions," |arXiv: 1102.2357 [hep-ex] 



ATLAS Collaboration, "Search for squarks and gluinos using final states with jets and 
missing transverse momentum with the ATLAS detector in ^/s = 7 TeV proton-proton 



collisions," |arXiv: 1102 . 5290 [hep-ex] 



CMS Collaboration, "Search for supersymmetry in pp collisions at ^/s = 7 TeV in 



events with two photons and missing transverse energy," arXiv: 1103.0953 [hep-ex] 



CMS Collaboration, "Search for physics beyond the Standard Model in opposite-sign 



dilepton events at = 7 TeV," [arXiv : 1 103 . 1348 [hep-ex] 



ATLAS Collaboration, "Search for stable hadronising squarks and gluinos with the 



ATLAS experiment at the LHC," arXiv: 1103. 1984 [hep-ex] 



J. Alwall, M.-P. Le, M. Lisanti, and J. G. Wacker, "Searching for directly decaying 



gluinos at the Tevatron," \Phys. Lett. B666 (2008) 34| |arXiv:080 3 . 0019 [hep-ph]] 
T. Sjostrand, S. Mrenna, and P. Skands, "Pythia 6.4 Physics and Manual," 



JHEP05 (2006) 026, arXiv:hep-ph/0603175, 



T. Sjostrand, S. Mrenna, and P. Skands, "A Brief Introduction to Pythia 8.1," 
\Comput . Phys. Commun. 178 (2008) 852^ ,arXiv: 0710. 3820 [hep-ph], 



See also http://home.thep.lu.se/~torbjorn/Pythia.html 



G. P. Salam and G. Soyez, "A practical seedless infrared-safe cone jet algorithm. 



JHEP 05 (2007) 086 


,arXiv: 0704. 0292 [hep-ph] 




See also http://projects.hepforge.org/siscone/ 



ATLAS Collaboration, G. Aad et al, "Performance of the ATLAS 6-tagging 
algorithms," ATL-PHYS-PUB-2009-018 (2009). 

CMS Collaboration, "Algorithms for b jet identification in CMS," CMS PAS 
BTV-09-001 (2009). 

CMS Collaboration, "Search for Randall-Sundrum gravitons decaying into two photons 
in 7 TeV pp colhsions with the CMS detector," CMS PAS EXO-10-019 (2011). 

D. A. Dicus, C. D. McMullen, and S. Nandi, "Collider implications of Kaluza-Klein 



excitations of the gluons," Phys. Rev. D65 (2002) 076007 |arXiv : hep-ph/OO 12259 



42 



